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O ■ Stochastic field equations for linearized gravity are presented. The theory is compared with the 

usual quantum field theory and qüestions of Lorentz covariance are discussed. The classical radia- 



tion approximation is also presented. 
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^ 1 Introduction 

Quantum systems which are described by a Schròdinger equation allow a stochastic inter- 
pretation. The Fenyes-Nelson model and its generalization {§-0 provide such a 

stochastic description of quantum mechanics. Several field theories have been considered in 
this context ]7fl— and the results have so far been of some interest. In the present paper, 
the generalized Fenyes-Nelson model is applied to the weak field approximation of Einstein's 
general theory of relativity, the so-called linearized gravitational field [ 12 1 — [ 14 . 

It is worthwhile to attempt a stochastic model of quantum gravity for several reasons. 
First, since the gravitational field in the classical theory has the interpretation of a mètric 
tensor, it is disturbing that in the usual quantum formulation this geomètric interpretation 
is completely obliterated because the field and its derivatives become abstract operators 
on a rigged Hilbert space. Because of this difficulty it has become fashionable to play 
down the geomètric role of gravity when dealing with quantization, as in the formulation of 



Weinberg WM. A stochastic formulation of the quantized gravitational field may revitalize 



the geometrical interpretation of quantum gravity. 

Second, the stochastic formulation of quantum mechanics, at least in spirit, is a progeny 
of Einstein's profound discomfiture with the complement arity vision of Bohr. Since Einstein's 
views were based at least in part on considerations of the general theory of relativity, it seems 
fitting to pursue a stochastic model of quantum gravity. 

Third, the great successes which gauge theories have had in the theory of elementary 
partides suggest that fundamental efforts such as the stochastic reformulation of quantum 
mechanics be concentrated in this general area. 
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The linear theory of gravity has been selected for analysis because it avoids the extremely 
difficult problem of divergences in the full theory, and it is sufficiently simple to allow, 
perhaps, the beginning of a probabilistic geomètric interpretation of quantum gravity. 

In Sec. II the linearized Einstein field equations are briefly recounted. In Sec. III the 
usual quantization procedure is outlined. In IV the method of stochastic quantization is 
applied to the linear theory, and in V the curious random classical radiation approximation 
is presented. 



2 The linear field equations 

The Einstein field equations in the linear approximation may be written 

where the notation of Ref. 13 is followed except that the gravitational constant G is not 
set to unity. Indices are raised and lowered via the Minkowski mètric n^. The full mètric 
tensor is related to the above fields by 

9fj,y = + h^, (2) 

h^ u 2^V^q' (^) 

h/xu = — —n l·lI/ h a , (4) 

where is the energy-momentum tensor of all nongravitational sources. 

Once the linear theory is solved, classically, the nonlinear corrections of the full theory 
may be included as a perturbation by introducing a suitable term in the energy-momentum 
tensor (Ref. 14, p. 165). This procedure has not worked so far in the quantum theory because 
of the celebrated divergences introduced by the nonlinear terms. 

Weinberg has suggested that general relativity be treated as an ordinary field theory with 
the linear theory as the starting point. This will be the approach taken here, although one 
of the primary motives for seeking a stochastic model of gravitation is to restore a geomètric 
interpretation to the quantum theory. 

The linear theory possesses a gauge invariance similar to electromagnetism (Refs. 13, p. 
439; 14, p. 254). The field equations are left invariant under the transformation 

V ~> V - dCJdx" - dC„/dx», (5) 

where C is an arbitrary 4-vector field. All observables are independent of the gauge, and 
so gauge conditions may be imposed by fiat in order to facilitate solution. We shall work in 
the Lorentz gauge where one requires 

K, a = 0. (6) 

The field equations become in this gauge 

- d a d a V = 167rGT M „. (7) 
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The solutions to Eq. (0) may be written as a sum of a retarded solution plus a free field 
solution, 



h — h R 4- h' m 



where 



and 



h* = / 4GT^(x',í-|x-x'|)/|x-x'|rfV (9) 



d a d a h™ = 0. (10) 

The gauge condition [Eq. @] does not determine the gauge completely. Any additional 
transformation such that 

d a d a C, = 0, (11) 

will preserve the Lorentz condition. This gauge transformation of the second kind can be 
used to make the free part of the field equation both transverse and traceless (Ref. 13, p. 
946). No generality is lost, therefore, by requiring 

Ku = 0, (12) 
K = 0. (13) 

Periòdic boundary conditions in the spatial coordinates will be imposed as an aid to 
quantization. If the length associated with this periodicity is L, then wave numbers in 
a Fourier decomposition are restricted to have components which are integral múltiples of 
2tt/L. Although the energy momentum tensor must also have this periodicity for consistency, 
no generality is really lost since L will be taken to infinity eventually. With these conditions, 
the Fourier decomposition of the free part of the gravitational field is 

KTj = (l/v / 2L 3 / 2 )^e ?J (A,k)e ikx Q(A,k,í), (14) 

A,k 

where we have dropped the distinction between h m and h m , since in the transverse traceless 
gauge they are the same. A denotes the polarization state and it can take on one of two 
vàlues. The e's must satisfy the conditions 

&ij ^ji} ki^-ij 0, SijSij 1, €.{{ 0, 

e^A) e^A') = Í A)V , e^(A, -k) = e*. (A, k). (15) 

The polarization tensors can be chosen real, and, with k in the z direction, they can be taken 
as 

= 4? I - 1 I > (16) 
€,-,■( x) — i u u | . (17) 
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for the two independent states of polarization denoted by + and x . The following expression 
is often useful: 

E ey(A)e„(A) = \ (-6*5* + 5*5* + 5*5%) , (18) 

where 

5* = 5 i3 - hkj/k 2 . (19) 
The equations of motion for the field require that the Q's satisfy 

Q + k 2 Q = 0, (20) 

so that a Lagrangian may be chosen of the form 

L = fI t (\\Q\ 2 - 1 ^ 2 \Q\ 2 ), (21) 

where / is an as yet undetermined parameter. In the sum over k each Q appears twice 
because reality requires: 

C(A,-k) = G'(A,k). (22) 

The conjugate mo menta are 

/ReQReP, -^U- = /ImQ = ImP, (23) 



OReQ dlmQ 
and the Hamiltonian may be written as 

H = E(V| 2 /4/ + Ík 2 /|Q| 2 ) 

A,k V 4 7 

= fJ2(\\Q\ 2 + ^ 2 \Q\ 2 )- (24) 

This must now be compared to the actual energy of the gravitational field to set the scale 
for quantization. 

The physical energy to be associated with a free gravitational field has been calculated 
(Ref. 13, p. 955). It is given by the 00 component of the stress energy tensor: 

^ = 3^gX(^)(^)^' (25) 
where the integration is taken over a three-dimensional cube (íl) of edge length L. One finds 

T 00 = (1/327TG) / {Kjfd^x 

= (l/3 27 rG)E Q|Q| 2 + Ík 2 |Q| 2 ), (26) 

so that we may make the identification 

/ = 1/327TG, (27) 

and therefore, 

P = (1/32ttG) Q. (28) 
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3 Quantization 

Quantization is achieved by imposing commutation rules as follows at equal times: 

(1/32ttG) [ReQ(A, k), ReQ(V, k')] 

= -ï^AA'(^k,k' + 5 k ,_k')> 

(1/32ttG) [ImQ(A, k),ImQ(A', k)] 

In terms of the fields, the equal time commutation rules become 

hij (x) , h k i (y) = -ih5 ijH (x - y ) , 



where 



and where 



8ijki(x - y) 



d 3 k 
(2^)3 



e ik(x - y) ^ fcZ (k), 



SijkiQt) = eii(A,k)e fcí (A,k) 

A 



(29) 
(30) 

(31) 
(32) 
(33) 



is given by Eq. (TE). The dynamical relation which fixes the quantum theory is 

Q = -k 2 Q. (34) 



The field theory is equivalent to uncoupled harmònic oscillators. In the ground state one 
finds 

d 4 k S--u,Ck) p iklí ( x »-y>*) 



(35) 



(27r) 4 -k^k^ + ie 
Using the rule of Gaussian combinatòries, 

(0 \T4> X {X X ) X ... X (f> n { Xn )\ 0) = £(0 \T<f>l{Xl)Mx2)\ 0) X ... X (0 |T0„-lK-l)0n(z„)| 0), 

(36) 

where n is even, T denotes time ordering, ir denotes a sum over distinct permutations, and 
the higher oder moments for the ground state of the quantum field may be calculated. The 
spectral energy density for the field is found to be the same as for electromagnetism: 



p(iü) = (h/2lï 2 )LÜ 3 



(37) 



corresponding to an energy of \ fiui for each degree of freedom of the filed in the ground state. 



4 Stochastic quantization 

The linearized gravitational field as a quantum system may be reduced to an infinite dimen- 
sional Schròdinger equation of the form 



ip = ih 



(38) 
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The dimensionality of this equation can be made finite by imposing a cutoff in momentum 
space. Since much of the theory upon which stochastic quantization relies has been done for 
only finite dimensional systems, it is convenient to impose a momentum cutoff which may 
be taken to infinity in most expressions of interest. This is not an important limitation and 
it can be expected as the theory of Markov fields advances it will eventually be eliminated. 
For the present, we shall assume a momentum cutoff. Let us chose the following notation 



E 



d 2 



dQÍ 



2 ; 



(39) 



where R and S are real functions. 

Direct computation shows that the following equation is equivalent to Eq. so long 
as R and S have a first time derivative and a second Q derivative, and so long as z has a 
nonzero real part 



(zkf 



A Q + ÍV(Q) 



2 



{* - 1) 



R N 



D R+iS/z 



i(zh) 



d_ 
dt 



R+iS/z 



(40) 



where z is a constant which may be complex. 
Suppose that z is is purely imaginary, so that 



%\z\. 



Then Eq. (|4"0|) is still true if (|38| ) is true, but (pEPP is only one real equation 



(\z\hf 



Ao+ f r-i- — (i + 



2 




R+S/\z\ 



\z\h 



d_ 

dt 



,R+S/\z\ 



(41) 



(42) 



This equation must be supplemented by another real equation since the original Schrodinger 
equation contains two real equations. Another equation may be generated by choosing 



-i\z\ 



which leads to 

\\z\hf 



( h 1 
a q + (v + — {i + \„ 



A Q e R * 



,RS/\z\ 



\ z \hL· a - 3 M. 

dt 



(43) 



(44) 



V 2 v ■ ■ / e 
taken together are equivalent to Schròdinger's equation by direct 



Equations (|42 ) and 
computation. 

Consider the two real Eqs. (f42|) and (0). They have the mathematical form of the 
heat equation. Equations of this genèric form are characteristic of a certain type of diffusion 
problem called generalized Brownian motion or Ito processes. In Ref. 4 equations of this form 
were derived within the context of diffusion theory. The fact that Schròdinger's equation 
can be rewritten in the form (^) and (|44] ) is the basis for the stochastic interpretation of 
quantum mechanics. 
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The stochastic interpretation rests on the hypothesis that all of the experimentally veri- 
fiable predictions of quantum mechanics may be deduced from Schròdinger's equation. This 
hypothesis shall be assumed true in this paper. 

Stochastic quantization is achieved, following Nelson |2|, || by associating with each co- 
ordinate Q a stochastic process which is defined by the stochastic differential equation: 

dQi = 6 f (Q J í) + £W i (í), (45) 

where the W's are Wiener processes which satisfy 

E{dWidWj) = 2v5 ij dt. (46) 



Processes defined by Eq. fl4"5| ) go by the name: generalized Brownian motion, Ito processes, 
or multidimensional diffusion processes. The mathematical background contained in Ref. 2 
is sufficient for an understanding of stochastic quantization. Other good sources for this 
subject are Refs. 15-22. Since the formalism of forward and backward derivatives has been 
developed only by Nelson D , a careful reading of his book is essential to an understanding 
of stochastic quantization. A method relying on an operator formalism, rather than the 
forward and backward time derivatives was presented in [|J, but it is completely equivalent 
to the Nelson procedure. 

Many existence and uniqueness theorems for the processes defined by Eq. (]4q) are pre- 
sented in Refs. 2, 15-22. The most important theorems show that if b satisfies a global 
Lipschitz condition then Eq. ( f45"D has a solution which is a continuous Markov process and 
which is unique (for example, Ref. 2, page 43). 

In Ref. 2 forward and backward time derivatives D and D* are defined by 

Df(Q,t) = lim l E (f(Q(t + h),t + h) 

-f(Q(t),t)\Q(t) = Q), (47) 

DJ(Q,t) = lim \E{f{Q{t),t) 

-f(Q(t-h),t-h)\Q(t) = Q), (48) 

and independent of any dynamical assumption, these satisfy 

DQi = k, (49) 
D*Qi = b u , (50) 

bi-K = 2uJ^ln\p(Q,t)}, (51) 
D = | + E^^A S , (52) 

D. = l + Ç^-^Aa. (53) 
The dynamical assumption which leads to Schròdinger's equation is 

\{DD* + D*D)Qi + f (D - D*) 2 Qi = -J^V, (54) 
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where (3 is a constant and where 



v = h/2y/(í-P/2). (55) 
The Schròdinger wave function is written in the form 

iP = e R+izSN , (56) 

where 

z = l/(l-/3/2) 1/2 (57) 

and 

h = 2v-^-{R + S N ). (58) 

It is straightforward to show that Eq. (153) implies that ( 56| ) satisfies (BSp. 
If z/ is to be real, then we must demand the condition 

(3<2. (59) 

It is possible, because of Eq. ([55]), to choose any value of the diffusion parameter v for a 
stochastic model of quantum mechanics. In the limit when v — > 0, the model becomes 
deterministic and equivalent to Bohm's hidden variable thoery ||23|| , as was first pointed out 
in i. 

In applying this formalism to linear gravity, we use the Q's in Eq. fl!H]), but mindful of 
the condition (|22|) which relates Q's for antiparallel wave vectors. If we define 



K = Dh§, (60) 
hj, = D.tifi, (61) 

then with the aid of the operator 



A,k ' d . d 



A iJ '(x, t) = Y e ij (X, k)e ïk ' x x — — + i — — , (62) 

{V2)L 3 / 2 ^ K ' \dReQ dlmQ j v ; 

the forward and backward time derivatives can be expressed as 

d 



D = —+ í b^Sx + v [ A ij A tí d 3 x, (63) 
ot Jn Jn 

D, = — — h í b ir A ij d 3 x~u í A ij A ij d 3 x. (64) 
Ot Jn Jn 



If we define a random field by 

WH^t) = ^ L3/2 J2e^(X,k)e^W XM (t), (65) 

where the real and imaginary parts of Wx^ are Wiener processes which are independent of 
one another and which satisfy 

E[d(ReWx, k )d(ReWx', k >)} = 2v6xx>(5u,v + S^. v )dt, (66) 
EidilmWxïWmWx',*)} = 2v6xx>(.&i*,u - &i*,-u)dt, (67) 
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then the stochastic differential equation may be written as 



dh% = bijdt + dWij. 



(6É 



With the gravitational field expressed as in Eq. (|8|), and with the retarded solution still given 
by Eq. (PI), the field equations for the free part of the field become 



~(DD* + D*D) + ]-/3(D-D*) 2 
Z o 



d,d l 



h 1IL 



0. 



(69) 



which is the stochastic version of Eq. ( |T0| ) in the transverse traceless gauge. 

In order to solve these equations, Schròdinger's equation must first be solved in Q space. 
Then the ò's are calculated using Eqs. fl56|) and (|58D. Once the ò's are known, Eq. fl68|) can 
in principle be solved. 

In general, it is more difficult to solve the stochastic equations than to solve the usual 
quantum mechanical equations. Even for the free field the stochastic processes for excited 
states are difficult to calculate. So far, the stochastic method has proven useful in practical 
problems only for stationary state problems where considerable simplifications occur. See, 
for example, Simón |24| for a review of results in this area. Although Simón does not 
explicitly make the connection with Nelson's theory, many of the methods he discusses 
may be considered as applications of the Fenyes-Nelson model to stationary state quantum 
systems. 

We now illustrate the theory for the ground state field. The solution to Schròdinger's 
equation in Q space for the ground state is 



VKQ) = II ex P 

A.k 



|Q(A,k)p 



íú 



hl28nG 



íú 



up to a normalization constant. The ò's are found from Eqs. ([56]) and ( [58D with Sn 
One finds for the ò's of Eq. (Bl) 



A,k 



bv(x,t) = -2z//(v / 2)L 3 / 2 x ^e ÏJ (A,k) e 41 " 

The stochastic differential equation becomes 

dh ij = b íj dt + dW íj . 
Using the property of the in Eq. ( |66"1) . 



íú 



mino 



fl(A,M). 



E[dW v (x,t x )dW k \y,t y )\ = 2vdtó ílkl (x-y 



(70) 
= 0. 

(71) 

(72) 
(73) 



where the delta function is that of Eq. fl32|) , the stochastic equations may be integrated to 
yield 



E(h l \x)h kl (y)) = h(16TrG) 



d 3 k 
(2tt) 3 



3 ik-(x-y) 



x e -(v/16iTGh)\\i\\t x -ty\ gijkln^ 



(74) 
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Since b 13 of Eq. (|7Ï|) is linear in the Q's, the process turns out to be Gaussian with zero mean 
so the covariance (|74]) determines all higher moments from the rule ([361). The covariance 
may also be written four-dimensional integral 

E(h ij {x)h kl {y)) = (327rG^) / — -- " K *{ e ik ^- y) x e ^/{m^G)\t x -t y \ _ ^ 

J (27r) 4 Kq + k 2 



In expressions (f73|)-(|75|) the infinite volume limit has been taken. Nelson's value for the 
diffusion parameter is 

v N = 16nGh. (76) 

Comparing Eq. (|75|) with (|35|), it may be shown that the stochastic covariance in (|75|) may 
be obtained from the quantum covariance in (|3~5"D by analytically continuing 



t x -> -i[v/(hl6nG)]t x , t y -> z[i//(ftl67rG)]í w . (77) 

This procedure of analytic continuation yields the Schwinger function. Using the rule fl36"l) , we 
obtain the general result: The moments of the stochastic theory are equal to the Schwinger 
functions of the quantum theory with the times scaled by the factor v(TilQirG). When 
Nelson's value for the diffusion parameter is chosen [Eq. (76)] the stochastic covariances 
become equal to the Schwinger functions. This result is similar to the results obtained in 
scalar field theory [7, 9, 10] and in electromagnetism [8, 9, 11]. 

Examining Eq. ([75]), we see that the covariance is not manifestly Lorentz covariant, even 
discounting for the fact that we have chosen a noncovariant gauge. Lorentz covariance is 
violated by more than just a harmless gauge transformation in Eq. (|75|). This is a surprising 
and perhaps paradoxical result which has been known for some time 0-['lï]. Despite this 



lack of manifest Lorentz covariance in the ground state, there is good reason to believe that 
the experimental predictions of the stochastic theory are consistent with special relativity 
and are in fact the same as ordinary quantum theory. The argument is as follows. Since 
for any real value of v we have a stochastic model for a given solution to Schròdinger's 
equation, and since we believe that this equation contains all of the experimentally verifiable 
predictions of quantum mechanics, then there is reason to think that it is impossible to 
measure the diffusion parameter. If this is true, then all of the experimentally measurable 
predictions of the theory (scattering cross sections, line spectra, etc.) must be constant when 
considered as analytic functions of v. Because of this, we may consider analytic continuations 
to complex v without affecting measurable predictions of the theory. When the stochastic 
theory is continued to 

v = ih(lQirG), (78) 
then one finds that the moments of the stochastic theory become the Green's functions of 



quantum field theory. This result was also found for the scalar field [|TÜJ and the electro- 
magnètic field The stochastic theory can be expressed in a mathematical form which 
is identical to ordinary quantum theory when v is continued to the value in Eq. ([T|) (or its 
complex conjugate). See, for example, Ref. 25 for a derivation of the operator formalism of 



quantum mechanics within this framework. The methods of J25| generalized easily to the 
present theory in Q space. 
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Stochastic quantization in the ground state for a field theory leads to moments which 
are analytical continuations to imaginary times of quantum moments. The stochastic in- 
terpretation suggests that a certain reality be attributed to the theory so continued. It is 
interesting that imaginary time continuations have played a surprisingly important role in 
modern physics. Euclidean field theory |2B| has led to advances in constructive field theory. 
Complex manifold techniques [27] have led to a deeper understanding of gauge theories. 
Bound state problems like the Bethe-Salpeter equations are often best solved by making a 
"Wick rotation" to imaginary times. Perhaps the results of stochastic quantization provide 
an explanation for the usefulness of imaginary time continuations. 



5 The approximation of random classical radiation 

In electromagnetism it has proven interesting in several applications to approximate the 
ground state of the quantum field by a superposition of classical plane waves with radom 
phases. This approximation has become known as random electrodynamics [^] and it has 
been found that the diffusion of charged partides in harmònic oscillator potentials and 
exposed to such radiation leads convincingly to Schròdinger's equation with the correct 
nonrelativistic Lamb shift ||29|| . We present this random phase approximation for the linear 
gravity theory in the hope that it may find similar uses and also for comparison with the 
stochastic theory. 

We first write the mètric perturbations as general solutions to the free field equations. 
We consider only the free field in the transverse traceless gauge 

h ij (*,t) = ^(À,k)e* kx Q(À,k ; t), (79) 

Q(A,k,í) = 32nGh/\k\x {cos {tot + 9^, k)]+ i cos {cot + 9 2 (\,k)]}. (80) 
Reality demands 

e 1 (\, -k) = e 1 {\, k), e 2 {\, -k) = e 2 {\, k) + n. (si) 

The #'s are all independent of one another except for the conditions in Eq. fl8"T|). They are 
random phases which take on vàlues fro to 2tt. The averaging process is carried out by 
integrating over these phases. One finds for the covariance 

E e (h ij (x)h kl (y)) = 16nGh J cos (u (t x - t y )) ^«(k)/|k|. (82) 

It is easy to show that this is equal to the symmetrized quantum expectation 

E e (hV(x)h kl (y)) = (0| Sym(/i« (x)h kl (y))\0), (83) 
where Sym denotes the symmetrization operation 

Sym(0 1 (s 1 ) x ... x <p n (x n )) = Y,—; M x i) x . . . x (j> n (x n ), (84) 

p n\ 

and where P denotes a sum over all permutations of the arguments of the fields. 
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It can be shown that in the infinite volume limit the higher the moments of the random 
phase average satisfy the Gaussian combinatoric rule [Eq. ([36])]. This is true despite the 
fact that the individual Fourier coefficients in Eq. ( |79|) are not distributed normally. The 
reason is the central limit theorem, as the field is a sum of an infinite number of independent 
random variables. Since the symmetrized quantum expectations also satisfy the Gaussian 
combinatoric rule, it follows that 

E e {h hh { Xl ) x ... x l·l·l n {x n )) = {0\Sjm(h iljl ( Xl ) x ... x h injn (x n ))\0) . (85) 

A more detailed derivation of this result has been given for electromagnetism by Boyer f2~8 | 



whose analysis can be applied with little modification to the linear gravity theory to derive 
Eq. (|85|) . It is important to realize that ( j85|) is not true for a one-dimensional oscillator or 
even for a finite dimensional oscillator. It can only be derived in the present case in the 
infinite volume limit. Thus the random phase approximation does not give a very detailed 
model of the quantum field accurate down to the level of a few normal modes. It is unlikely, 
in the author's opinion, that a consistent interpretation of the full quantum field theory 
could be obtained from a classical random phase approximation for field theory. This is in 
sharp distinction to the stochastic quantization theory of Sec. IV which allows a consistent 
reinterpretation of all quantum phenomena. Still, the random phase model can be useful 
when considering the effects of vacuum fluctuations on matter. 
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